Complex Dynamics of Piecewise Linear Model of the Van der Pol Oscillator 

Under External Periodic Force 
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The electronic model of Van der Pol oscillator with piecewise linear V — /-characteristics of nonlinear 
element is proposed. It is carried out the experimental investigations of behaviour of the Van der Pol 
oscillator under external periodic force. The parameter plane "the amplitude of external force - the 
frequency of external force" for the dynamical regimes of the oscillator is experimentally plotted. 
It is shown in experiments and numerical simulations that transition to chaos in the oscillator 
takes place through the period doubling cascade. It is revealed that the internal structure of the 
synchronization tongues has the characteristic view as for the crossroad area. 

PACS number(s): 05.45. +b 
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The Van der Pol oscillator is a classical example of 
self-oscillatory system. More than seventy years after 
its appearing, the oscillator has become the main model 
self-oscillatory system with two dimensional phase space 
[1-3]. The classical experimental setup of the system is 
the oscillator with vacuum triode [2,3]. The mathemati- 
cal model for the system is well known second order or- 
dinary differential equation with cubic nonlinearity - the 
Van der Pol equation. The Van der Pol equation is ba- 
sic model for oscillatory processes in physics, electronics, 
biology, neurology, sociology and economics [4]. 

It was much attention dedicated to investigations of 
the peculiarities of the Van der Pol oscillator behaviour 
under external periodic (sinusoidal) force and, in partic- 
ular, the synchronization phenomena and the dynamical 
chaos appearing (see e.g., [5-9]). 

The investigations of the forced Van der Pol oscilla- 
tor behaviour have carried out by many researchers. So, 
in 1927 Van der Pol and Van der Mark [1] in their in- 
vestigations of the oscillator behaviour in the relaxation 
oscillation regime found that the subharmonical oscilla- 
tions are appeared during changes of natural frequency 
of the system. Moreover, the authors noted appearing 
of "irregular noise" before transition from one subhar- 
monical regime to another. It seems to be it was one of 
the first observations of chaotic oscillations in the elec- 
tronic circuit. Later, in 1945, Cartwright and Litlewood 
[10] during their analyzing the Van der Pol equation with 
large value of nonlinearity parameter have shown that the 
singular solutions exist. In 1949 Levinson [11], analyti- 
cally analyzing the Van der Pol equation, substitutes the 
cubic nonlinearity for piecewise linear version and shows 
that the equation has singular solutions in the case also. 
It was revealed that chaotic behaviour is appeared in the 
Van der Pol equation with smooth nonlinearity by period 
doubling bifurcations [9] and crossroad area is observed 
in synchronization tongues [12]. 

However, up to the present day the difficulties con- 



nected with creation of the electronic model of the nonau- 
tonomous Van der Pol oscillator which demonstrates 
transition to chaos in physical experiment and will al- 
lows adequate numerical simulation are remained. The 
efficient enough way for creation of such kind models is 
building of the systems with piecewise linear voltage- 
current characteristics [13,14]. In the paper [14] on the 
basis of numerical simulation shown that in piecewise lin- 
ear model of sinusoidally forced Van der Pol oscillator 
the transition from periodic oscillations to chaotic ones 
is possible. The transition takes place through cascade 
of period doublings. 

In this paper the results of experiment and numeri- 
cal simulation of the sinosoidally forced Van der Pol os- 
cillator model with symmetric piecewise linear voltage- 
current characteristic of nonlinear element are discussed. 

Experimental circuit and models. The classical dimen- 
sionless form the nonautonomous Van der Pol oscillator 
is described by equation 



e(l — x 2 )x + x = acoswt, 



(1) 



where x is the state variable, e is the nonlinearity param- 
eter, a and to are the amplitude and frequency of external 
force, respectively. 

For building of electronic circuit of the Van der Pol 
oscillator which would be adequate for modeling by (1) 
the serial circuit could be used (Fig. la). The nonlin- 
ear element (— r) in the circuit (Fig. la) must obey with 
symmetrical characteristics with negative resistance re- 
gion. The nonlinear resistor with three-segment piece- 
wise linear S 1 — type voltage-current characteristic is cre- 
ated on an operation amplifier [15]. The scheme is shown 
in Fig. lb, where the rest parameters are Ri=R3=3.6 kfi, 
R 2 =400 n, E+ = 12 V, E_ = -12 V. The voltage- 
current characteristic measured in experiment for the 
above mentioned parameter values is presented in Fig. 2. 
The resistance at negative region of nonlinear element 
voltage-current characteristic is R- = —400 fl and at 
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positive region is R + = 3600 £1. 

On the basis of Kirhgoff laws the equations for the 
Van der Pol oscillator model (Fig. la) could be written 
as 



L 



dl 



dt 
dV c 
dt 



V c + RI + V(I) = A cos cot, 
I 

C' 



(2) 



where Vl is the voltage at the inductor, Vr is the volt- 
age at the resistor, Vc is the voltage at the capacitor, 
V(I) is the function describing piecewise linear voltage- 
current characteristic of the nonlinear clement " — r" , 
A and to are the amplitude and frequency of external 
force, respectively. Introducing new dimensionless quan- 
titles: r = co t, x = I/Jo, y = V/V , f(x) = V(I)/V , 

= {V /I )y/UjL, A = 



fl = Lo/cOq, UJq = l/\ 



(A/I ) y/CjL, G = R^/Cjl in term of which (2) be- 



come 



x = —Gx — ey — ef(x) + A cos fir, 
1 

y = -x, 



(3) 



where Jo and Vb are the values of current and voltage in 
extremum of the voltage-current characteristic (Fig. 2), 
(•) = d/dr. 
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FIG. 1. Schematic diagram of the Van der Pol oscillator 
under external sinusoidal force (a) and the schematic diagram 
of nonlinear element with piecewise-linear characteristic (b). 

Along with the piecewise linear form of f(x) in 
the equations (3) the arbitrary function f(x) could be 
used. So, if nonlinearity in eqs. (3) is governed by 
f(x) = ^x 3 — x then additive " —Gx" could be removed 
by rescaling and finally the eq. (1) is obtained. Thus, the 
experimental scheme (Fig. la) and the appropriate equa- 
tion (3) are adequately described the classical Van der Pol 
oscillator under external force (1). 
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FIG. 2. The piecewise-linear voltage-current characteristic 
of nonlinear element which depicted in Fig. lb. 

Experimental results. For systems under external force, 
as a rule, the synchronization phenomena of the natural 
oscillations by external periodical force are investigated 
at the parameter plane "the amplitude - the frequency 
of external force" . Most frequently the synchronization 
regimes are marked by (1:1), (1:2), (1:3) and so on. Such 
notations of the synchronization regimes (m : n) show 
the relation of the natural frequency of gerenator's os- 
cillations (m) to the frequency of the external force (n). 
There are so called main synchronization regimes marked 
by (1:1), (1:2), (1:3), (3:1), (2:1) etc. Between such 
regimes the more complicated multiloop regimes of high 
order are placed. All these regimes have the rational 
relation of frequencies. Moreover, the regions (quasiperi- 
odical oscillations) with irrational relation of frequencies 
exist. 

Experimental investigation of the driven Van der Pol 
oscillator behaviour is in plotting of synchronization 
regimes at the parameter plane. For experiments in this 
work the following parameters were used: i?=25.8 f2, 
C7=83.5 nF, L=4.37 mH (Fig. la). The resistance R 
is total resistance (excluding —r, Rl—23.3 ft is the re- 
sistance of inductor) of the circuit. The boundaries of 
oscillatory regimes were defined from the phase portrait 
view during the changes of controlling parameters. In the 
some cases the changing of both parameters simultane- 
ously is needed because the oscillator demonstrates the 
multistable behaviour. The multistability is appeared as 
existence of the histeresis boundaries for different oscilla- 
tory regimes. It could be proposed there are the regions 
which could be reached for the both parameter changing 
only. 

When the parameter plane has been plotted the syn- 
chronization regimes were identified by phase portraits. 
The different regimes depicted as the synchronization 
tongues. Generally speaking, the tongues of higher or- 
der are placed under the tongues of the lower order, for 
example, the synchronization tongue (1:2) could be par- 
tially placed under the synchronization tongue (1:1). The 
regions with period doubling is marked with subscript 
letters (e.g., (1 : 2) 2 , (1 : 2)4 ), the regions with chaotic 
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behaviour is marked with subscript C (e.g., (1:2) c)- 

The parameter plane "the amplitude - the frequency of 
external force" for the Van der Pol oscillator (Fig. 1) with 
the piecewise-linear ^-/-characteristics of nonlinearity 
(Fig. 2) is shown in Fig. 3. The main synchronization 
tongues are shown in Fig. 3 only, although, many narrow 
tongues with rational frequencies relation are observed 
in experiment. The phase portraits for typical dynami- 
cal regimes depicted at the parameter plane (Fig. 3) are 
shown in Fig. 4. In Fig. 4 the amplitude of external force 
changes for the horizontal axis and the voltage in the 
point of connection of capacitor C and inductor L - for 
vertical axis. The data were sampled with an analog to 
digital converter with 8-bit resolution at sampling rate 
2 • 10 7 samples per second, the length of time series is 
3000 points. 
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FIG. 3. The parameter plane "the amplitude of external 
force - the frequency of external force" (^4 = \/~2 « 1.414 A', 
uj = 27r/o) for the Van der Pol oscillator (Fig. 1) with piece- 
wise-linear ^-/-characteristic of nonlinear element (experi- 
ment) . 

As it is seen in Fig. 3 the overlapping of the syn- 
chronization tongues (1:2) and (1:3) in the vicinity /o ~ 
12.6 kHz, A' ~ 0.6 V takes place. The chaotic regimes 
with strange attractor could be seen in the synchroniza- 
tion tongues (1:2) and (1:4), more precisely in (1 : 2) 4 , 
(1 : 4) 4 regions. It is easily seen that the internal struc- 
ture of the synchronization tongue (1:4) has the view that 
is characteristic for the bifurcation composition called 
crossroad area [16]. 

Numerical simulation. In the numerical simulation of 
the driven Van der Pol oscillator with piecewise-linear 
characteristic the dimensionless form of equations (3) is 
used. The voltage-current characteristic of the nonlinear 
element (Fig. 2) is transformed to dimensionless form 
also (Fig. 5). The analytical form of piecewise-linear V- 
/-characteristic could be described by equation 



9.09a; + 10.09, x < -1.0, 
f(x)={-1.0x, -1.0 < x < 1.0, 

9.09a; - 10.09, x > 1. 



(4) 



Taking into account the values of the characteristic 
scales of the piecewise linear nonlinearity (Jo = 3.0 mA, 
Vb = 1.2 V) and the values of the circuit elements (Fig. 1) 
the dimensionless parameters (G = 0.113, e = 1.75) of 
the system (3) are calculated. 
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FIG. 4. The phase portraits for A'=0.60 V and 
fo = 6.5 kHz (1:1) (a), f = 11.7 kHz (1:2) (b), fo = 12.0 kHz 
(1 : 2)2 (c), f = 12.15 kHz (1 : 2) 4 (d), fo = 12.25 kHz (l:2) c 
(e), fo = 24.7 kHz (1:4) (f) for the characteristic dynamical 
regimes depicted at the parameter plane (Fig. 3). 

The parameter plane "the amplitude - the frequency 
of external force" for the system (3) with piecewise- 
linear characteristic (4) obtained in numerical simula- 
tion is shown in Fig. 6, where synchronization regions 
are marked in gray scale. For plotting of the parameter 
plane (Fig. 6) the system (3), (4) was integrated by using 
of the fourth order Runge-Kutta method with constant 
time step h — 0.05. The identification of the synchro- 
nization regime was carried out from Poincare section by 
selecting points in the section through the period of the 
external force. 

It is seen from Fig. 6, that the characteristic syn- 
chronization tongues are obtained. Main synchroniza- 
tion tongues are marked, besides, many synchronization 
tongues of higher order are also visible. The route to 
chaos via period doublings is observed within the syn- 
chronization tongues (1:2) and (1:4). Moreover, the view 
of the internal structure of the synchronization tongues 
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(1:2) and (1:4) allows to conclude that the crossroad area 
takes place in a similar way as it was found in natural 
experiments with the electronic circuit (Fig. 1). 

The visual comparison of the parameter planes ob- 
tained in the natural experiment and the numerical sim- 
ulation (Fig. 3 and 6, respectively) leads to conclusion 
about qualitative agreement of both experiment and sim- 
ulation. The further comparisons could be done if to 
rescale the experimental results (Fig. 3) to numerical sim- 
ulated ones (Fig. 6) using the coefficients introduced for 
obtaining the dimcnsionless equation (3) then it is seen 
the good quantitative agreement of both experiment and 
simulation. 
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FIG. 5. The rescaled piecewise-linear voltage-current char- 
acteristics of the nonlinear element (Fig. 2), described by 
eq.(4). 
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FIG. 6. The parameter plane "the amplitude of external 
force - the frequency of external force" for the Van der Pol 
oscillator with piecewise-linear V-J-characteristic of nonlin- 
ear element (numerical simulation). 

Conclusion. In this paper the electronic model of the 
Van der Pol oscillator with piecewise linear character- 
istic of nonlinear element is proposed. It is discussed 
results of the experiment on the investigation of the os- 
cillator behaviour under external sinusoidal force at the 
parameter plane "the amplitude - the frequency of ex- 
ternal force". The primary attention has been paid to 
the revealing of regions on the parameter plane in which 



the chaotic behaviour is realized and to investigations of 
internal structure of such regions. It is shown that tran- 
sition to chaos via period doubling is taken place within 
the synchronization tongues (1:2) and (1:4) (here ratio 
of natural frequency of oscillator to external force fre- 
quency) as well as in many synchronization tongues of 
higher order. 

The comparison of the parameter planes obtained in 
experiment and numerical simulation showed that the 
planes obey with both qualitative and quantitative sim- 
ilarity for boundaries of regions with synchronous and 
chaotic behaviour. To demonstrate this fact the typical 
points from the experimental parameter plane (Fig. 3) 
could be rescaled to dimensionless case and compared 
with the same points for the parameter plane obtained 
in numerical simulation (Fig. 6). So, for the point 
(f = 12.07 kHz, A' = 0.55 V) of minimum for A' at 
the boundary of (1:2) and (1:2) 2 regimes (Fig. 3) could 
be calculated dimcnsionless values (17 = 1.45, Ao = 1.13) 
that in good agreement with the point of minimum for Aq 
at the same boundary (Fig. 6), analogously, for the point 
(/o = 25.72 kHz, A' = 0.63 V) of minimum at the bound- 
ary (1:4) 2 and (1:4) 4 is obtained (17 = 3.08, A Q = 1.30). 
The region at the parameter plane where transition to 
chaos takes place is displaced within the synchronization 
tongues (1:2) and (1:4). The region has the crossroad 
area like view. The results of experimental investigations 
and numerical simulations of the changes of the parame- 
ter plane structure for large values of external force and 
the circuit elements will be published elsewhere. 
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